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A relativistic neutral scalar field is investigated in non-equilibrium thermo field dynam- 
ics. The canonical quantization is applied to the fields out of equilibrium. Because the 
thermal Bogoliubov transformation becomes time-dependent, the equations of motion for 
the ordinary unperturbed creation and annihilation operators are modified. This forces us 
ly^ | to introduce a thermal counter term in the interaction Hamiltonian which generates addi- 

tional radiative corrections. Imposing the self-consistency renormalization condition on the 
total radiative corrections, we obtain the quantum Boltzmann equation for the relativistic 
scalar field. 

Qh §1. Introduction 



A scattering process is described by a relativistic quantum field theory at high 
energy. The quantum field theory has to be extended to study an out of equilib- 
rium system for a relativistic field. There are several formalisms to introduce a 
thermal property in the quantum field theory, for example, thermo field dynamics, 
closed time path formalism and Langevin equation.^"® The thermo field dynamics 
(TFD) is a real time formalism based on the canonical quantization.®^ T n XFD the 
ly-^ ' thermal Bogoliubov transformation is introduced with particle number density and 

the thermal Bogoliubov transformed oscillator defines the so-called thermal vacuum 
state. The thermal average of a dynamical operator is represented as an expectation 
value in the thermal vacuum state. 

Much attention has been paid to generalize it for an out of equilibrium system 
^ \ and derive the time evolution equation for a particle number distribution. The 

non-equilibrium thermo field dynamics (NETFD) is constructed by extending the 
thermal Bogoliubov transformation in Ref. 8). The Boltzmann- like equation for a 
non-relativistic particle is found by the time evolution equation for the expectation 
value of the particle number operator with perturbed oscillators in Ref. [9]). 

In NETFD the thermal counter term is introduced from the consistency of the 
time evolution for an ordinary oscillator and the thermal Bogoliubov transformed 
one I3SHI2J) rj^g Boltzmann equation is derived from the self-consistency renormal- 
ization condition which is imposed on the quantum corrections with the thermal 
counter termP^'^ The method was extended to an inhomogeneous system with 
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diffusion IISiED Recently NETFD is applied to evaluate the Bose-Einstein conden- 
sation in trapped cold atom systems P^J"^ 

NETFD for a relativistic field is necessary to study the thermalization process 
of quarks and gluons in ultrarelativistic heavy ion collisions. An out of equilibrium 
system for a relativistic field is also essential in critical phenomena at early universe. 
Several works have been done for a relativistic field in NETFD.^D 1 '^ However, the 
canonical formalism of NETFD has not been fully established for a relativistic field 
yet. Thus we have launched our plan to make a systematic study of the canonical 
quantization in NETFD based on the thermal Bogoliubov transformation. 

In this paper we focus on a relativistic neutral scalar field and investigate the 
canonical quantization in a homogeneous system. In §2 we briefly review TFD for a 
scalar field. The thermal Bogoliubov transformation is made time-dependent for an 
out of equilibrium system in NETFD, in accordance with time-dependent particle 
number distribution. In §3 we discuss how to decompose the neutral scalar field. In 
§4 the scalar field is quantized in the canonical formalism. We introduce the thermal 
counter term and calculate the scalar propagator. The self-consistency condition is 
introduced from the structure of the scalar propagator. The time evolution equation 
is obtained from the self-consistency renormalization condition. In §5 it is confirmed 
that the self-consistency condition implies the correspondence between the thermal 
Bogoliubov parameter and the particle number density, given by an expectation of 
the Heisenberg number density operator. We impose the condition on the neutral 
scalar field with self- interactions. Concluding remarks are given in §6. 

§2. Non-equilibrium thermo field dynamics 

There are several real time formalisms to introduce thermal dynamics in a quan- 
tum field theory. TFD is one of the real time formalisms based on the canonical 
quantization. In TFD the statistical average of an observable quantity is repre- 
sented as an expectation value of the observable operator under a pure state called 
thermal vacuum. The ordinary oscillator operators, after each degree of freedom is 
doubled as below, are related to new oscillator ones through a Bogoliubov transfor- 
mation whose annihilation operator annihilates the thermal vacuum. The thermal 
Bogoliubov transformation is made time-dependent to deal with a non-equilibrium 
state in NETFD. In this section we briefly review the framework of NETFD. 

In TFD the thermal degree of freedom is introduced by doubling the Fock space 
which is spanned by ordinary annihilation, creation operators, (a a'), and tilde ones, 
(a a'). For a bosonic field the operators satisfy the commutation relations, 

k), (2-1) 
fc), (2-2) 
0. (2-3) 

We define the tilde conjugation rules by 



[a p ,al] = (2nfd^(p- 
[a p ,al] = 0, [a p ,d k ] = 



(AiA 2 )~ = AxA 2 , 



(24) 
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(ciAi + c 2 A 2 r = c\Ax + 4A 2 , (2-5) 
(A)~ = A, (2-6) 

where Ai and A 2 mean any operators, and c\ and c 2 are c-numbers. 

Time evolution for the tilde operator is generated by the tilde conjugate of a usual 
Hamiltonian. The tilde-Hamiltonian, H, is constructed only from the tilde-operators. 
The time evolution of whole a physical system is described by a hat-Hamiltonian, 

H = H — H. (2-7) 

To develop the canonical quantization in NETFD it is more convenient to rewrite 
the annihilation and creation operators by the transformation called the time-dependent 
thermal Bogoliubov transformation, 



${t) = B(n p {t)) a Pa${t), (2-8) 
^(t) = a^t)B-\n p ( 

where the upper indices, a and (3, represent the thermal doublets, 



m) = 4(t)B-\n p {t)f«, (2-9) 



a P = 1 a* ) ' a r> = \ a P ~®p ) ' ( 2 ' 10 ) 

£ = (!)> £ = (3 -&)■ (2 ' n) 

The thermal Bogoliubov matrices, B and B~ l , are chosen as 

BMO) = ( 1 + _7 (4) -"[ (4) ) , (2-12) 

"■'("-('))=(! (243) 

so that the Dyson- Wick formalism can be used.®' Since we work in the interaction 
picture, the operators, £ p (t) and a p (t), depend on time. In a particular case of 
equilibrium, the Bogoliubov parameter n p is time-independent and is taken to be 
the Bose-Einstein distribution. The time dependence of the transformed oscillators 
is given by 

(2-14) 

&{t x )=&e^ : (2-15) 



where uj p describes the on-shell energy for the scalar field, uj p = \/p 2 4- m 2 . Though 
the energy, u p , is not always time-independent, we confine ourselves to a system 
with a time- independent oj p , for simplicity. It should be noted that non-trivial time 
dependence is induced for the oscillators, a p and a p , through the time-dependent Bo- 
goliubov parameter with momentum index, n p (t). For a homogeneous and isotropic 
system which we assume in our practical calculations below, n p (t) is a function of 
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the time and the magnitude of the momentum. As will be seen in §5, the correspon- 
dence between n p (t) and the particle number density, given by an expectation of 
the Heisenberg number density operator, is shown by imposing the self-consistency 
renormalization condition for a non-relativistic system.® JHLO 
The thermal vacuum is defined by 

Z P \0) = £ P \0) = 0, (2-16) 

ml = ml = o. (2-i7) 

The thermal vacuum state is invariant under the tilde conjugation, = \0), and 

= {6\. After we perform the thermal Bogoliubov transformations (\2-8\i and 
(|2-9p on the oscillators, £ and £, Eqs. (|2-16p and ()2T7p are rewritten as 

(1 + n p (t))a p (t)\6) = n p (t)al(t)\6), (2-18) 

n p (t)al(t)\0) = (1 + n p (t))a p (t)\0), (2-19) 

(9\a p (t) = (0\al(t), (2-20) 

(0\aUt) = (9\d p (t). (2-21) 



Note that the thermal bra and ket vacua are not symmetric due to the choice of 
the thermal Bogoliubov matrices in Eqs. (|2-12p and (|2-13p . From Eqs. (|2-20p and 
()2-2ip we can show that the interaction hat-Hamiltonian, Hint, satisfies the following 
condition,® 

(0\H int = 0, (2-22) 

(but Hi n t\0) 7^ 0), which enables us to use the Dyson- Wick formalism. As is known, 
the Bogoliubov transformation keeps the commutation relations. Thus the trans- 
formed operators, £ and £, satisfy the commutation relations. 

[£ P 4] = (27r)V 3 )(p-fc), (2-23) 

[i p ,il] = (27rf5^(p-k), (2-24) 

fe>4] = 0, & s ffc] = 0. (2-25) 

A scalar field is quantized under these commutation relations. It should be noticed 
that physical observables are constructed by the original operators, a and a'. We 
would like to evaluate the expectation value for such operators under the thermal 
vacuum, \0). 



§3. Relativistic neutral scalar field in NETFD 



In the previous section we have introduced the bosonic operators, a p and £ p . 
A scalar field can be represented by either of these operators. Since the thermal 
Bogoliubov transformation does not depend on the time variable in equilibrium, the 
time dependence of both the operators is described by the same Hamiltonian in the 
interaction picture. On the other hand, the time-dependent Bogoliubov transfor- 
mation induces a discrepancy between the time evolution for a p and £ p in NETFD. 
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We introduce a decomposition of the scalar field in terms of the operator a p in a 
consistent manner. 

A neutral scalar field can be represented by the operators, £ p and £ p , in the 
interaction picture, 

^ {X) ~ J ^^^^^ + ^ P {t*) T Ye- ip - X l (3d) 

^ {X) = J (2^3 J^ p ^ e ~ iP ' X + (Ut,) T r 3 ) a e^ x ), (3-2) 

where t% is the third Pauli matrix acting on thermal indices. The canonical conju- 
gates for the fields, <f>^ and are given by 

^ {X) = H) j J^]^^(t x yP- x - (n^fre-^l (3-3) 

W = H) / ^s]f^(-^ e ' iP ' X + fe^f^V^). (3-4) 
These fields satisfy the canonical commutation relations at the equal time, 

[^(t,x),^(t,y)] = i5^(x-y)rf, (3-5) 

[^(t,x),^(t,y)] = i5^(x-y)5^. (3-6) 

It is straightforward to quantize the scalar fields, <j>£ and (f>g, under the thermal 
vacuum in a canonical formalism. However, a canonical formalism is necessary for a 
scalar field written by the operators a and a in order to evaluate physical observables. 

The time dependence of the operators, a p and a p , can be fixed by the time 
evolution equations for the operators, £ p and £ p . Differentiating Eqs. (|2-14|) and 
(|2-15p with respect to the time variable, we obtain the time evolution equations for 
the operators, £ p and £ p , 

dtX(t x ) = ~ico p $(t x ), (3-7) 
dtj^t x ) = iu> p $(t x ). (3-8) 

Applying the thermal Bogoliubov transformation (j2-8j) for the operator £ in Eq. (j3-7j) . 
we obtain a time evolution equation for the operator, a, 

dt x a^(t x ) = -i(u p - ih p (t x )T )^a^(t x ), (3-9) 

where the matrix, Tq, is 




(3-10) 



We note that T 2 = 0. From Eq. ([Ml) it is found that the energy eigenvalue for the 
operator, a, relies on the time derivative of the thermal Bogoliubov parameter and 
written as 

Qf{t x ) = UpS * - ih p {t x )Tf. (3-11) 
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Hence we write the time evolution for the positive frequency part of the scalar field 
as 

rt. 1 "P 



a p t (t x ) = expl-i J dt s n p {t s )\ og. (3-12) 



We obtain the differential equation for the operator, a, by the thermal Bogoliubov 
transformation (|2-9p for Eq. f|3- 8j) . 

dt x <ip(tx) = a§(t x )i(u> p - in p {t x )T f a . (3-13) 

The time dependence for the negative frequency part of the scalar field is represented 
as 

r i' tx ~\ 0a 

a$(t x ) = a$expUJ dt s Q p {t s )\ . (3-14) 

The positive and negative frequency parts rely on the same energy eigenvalue. 

Both the operators, a p and d p are organized into the positive and negative fre- 
quency parts of neutral scalar fields, 

/ Wf^ p W^ P ' X + (^sfr e~ ip - X ) , (3-15) 

= / {a a p (t*)e- ip - X + (a p (t x fr 3 r^ x } . (3-16) 

Due to the non-Hermiticity of Eqs. (|3T2p and (|3-14p . the neutral scalar fields (|3-15|) 
and (|3-16p are not invariant under the time-reversal transformation. The canonical 
conjugate fields, 7r" and 7f°, are decomposed into 

d 3 p , j p , x _ ( -j. A T,a p -ipx^ 



(x 



< (*) = H) / (i^s V 2 R(^) e** - (r 3 a P (^) T ) a e-**} , (3-17) 



Using the commutation relations (|2-3|) . we calculate the equal-time commutation 
relations for (p a and 7r a and get 

=i5^(x-y)rf, (3-19) 

[^(t,a;),7rf(t,y)] =W( 3 )(x-y)<^. (3-20) 

Thus we obtain a decomposition for the scalar field by the operators a and a with 
an ordinary canonical commutation relations. 

We construct the Hamiltonian, Hq, which describes the time evolution of the 
field, (p a . The equations of motion for the fields, <p a and ir a , are derived from 
Eqs. QIEI2D and pT4"]) . 

f l + i r ^A T \ d tm ^(x)=irt(x), (3-21) 
= - ( 1 - i^^To) (-V* + m 2 )^(x), (3-22) 
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where we make the definition, 6jy x = \/—V 2 , + m 2 . Thus the Hamiltonian, Hq, is 
found to be 



He 



d 3 x 



1 



l-^T ) a '(-V^ m ^f(,) 



(3-23) 



The equations of motion (|3-2ip and (|3-22p are reproduced from this Hamiltonian. 

§4. Self-consistency renormalization condition 

In the previous section the Hamiltonian, Hq, is obtained from the equation of 
motion for the neutral scalar field, <fi a . Here we regard the Hamiltonian, Hq, as the 
unperturbed part and adopt the perturbation theory. In the thermal doublet nota- 
tion the quantum field theory for a neutral scalar field is defined by the Hamiltonian, 



H — Hq + H, 



int ■ 



(4-1) 



where Hq and Hi n t represent the free and interaction part of the hat-Hamiltonian 
for neutral scalar field, respectively. The free hat-Hamiltonian is written as 



Hn 



d 3 x 



The unperturbed hat-Hamiltonian, Hq, is given from Eq. (|3-23p by 

Hq = Hq — Q, 

where Q is called the thermal counter term and is found to be 



(4-2) 



(4-3) 



Q 



d s x 



2 uj Vx 



2 uj Vx 



(4-4) 



So the interaction hat-Hamiltonian in NETFD, denoted by Hj, is not Hj n t but has 
to include the counter term, 



Hi — Hi n t + Q, 
and the total Hamiltonian is rewritten as 

H = Hq + Hi. 



(4-5) 



(4-6) 



Below we develop the perturbation theory in NETFD with respect to the unper- 
turbed Hamiltonian, Hq, and the interaction Hamiltonian, Hj. 
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According to the Dyson- Wick formalism, we evaluate a quantum correlation 
function. Thus the scalar propagator is given by 

Df{t x ,t y ,x-y) = {0\T[<%(x)%(y)u(oo,-oo)]\0), (4-7) 

where T denotes time-ordering and the operator u(t, t') is given by 



u(t, t') = exp 



dt s Hj(t s 



(4-8) 



Below we drop the momentum label, p, in the propagator and the scalar field for 
simplicity. From (|2-22p the thermal bra vacuum is invariant under the time evolution, 
(0\u(t,t>) = (9\. 

The propagator (|4-7p is rewritten in terms of the transformed operators, £ and 

I 

D^f(t x ,t y ,x- y) 

= B-\n\ Vx \(t x )) a ^ [^-^(el^^H^.g^Cj/H^.-cjo)^) 
+ 0(t y - t x )(e\4%_(y)u(t y , t x )4>l + (x)u(t x , -oo)\9)] B(n^Jt y ))^ 

+B- 1 (n lVxl (t x )) a ^ [9{t x - t y )(9\4>l + (x)u(t X: t y ){^ + (y)T 3 }^u(t y , -oo)\0) 
+ 9{t y - t x )(9\{^ t+ (y)r 3 y 2 u(t y , t x )</% + {x)u(t x , -oo)\9)] {B~ l (n^ y ft y )) T r 3 }^ 

+{r 3 B(n ]Vx{ (t x )) T } a ^ [9(t x - t y )(9\{n<t>^(x)}^ U (t x , ty)4%_ {y)u(ty,-oo)\0) 
+ 9(t y - t x )(9\^_(y)u(t y ,t x ){T 3 <P^(x)}^u(t x , -oo)|0>] B(n^Jt y ))^ 

+ {r 3 B(n lVxl (t x )f} a ^ [e(^-^)(fl|{7S0e,-(s)P«(tx,ty){4+(l/)73P"(*w.-oo)l^> 
+ 9(t y - t x )(9\{^ + (y)r 3 y 2 u(t y , t x ){T 3 ^(x)}^u(t x , -oo)\9) 
x{B-\n^Jt y )) T T 3 }^, (4-9) 

where the fields 4>£ ± and (j>£ ± show the positive and negative frequency parts, re- 
spectively, 



4>l + {x) 



(2vr) 3 ^/%~J~ p 
d 3 p 1 



{Er 3 ) a e -• ■< 



(2^)3^^ ^ p 
(2tt) 3 ^FET 



d 3 P 1 



(2tt) 3 



iupt x —ip-X 



(4-10) 
(4-11) 
(4-12) 
(4-13) 
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From Eqs. ()2- 16j) and (|2-17p it is found that the propagator (|4-9p has the following 
structure with respect to the thermal Bogoliubov matrices, B(n p (t x )) and B(n p (t y )), 

D a J i (t x ,t y , X -y) = B-\n lVxl (t x ))(<f y) ) ^Jt y )) 

+B hv,|(y)^i (X)y) )B (n^Jt y )) r 3 

^■1^ ( Sty) df(x,y) ) B-\n^ { t y ))^ (4-14) 
where 

dr 2 (x,y) = 9(t x - t y )(9\^4x)u(t x ,t y )^_(y)u(t y , -oc)\6) 

+ 9{t y - t x ){6\^lSy)u{t y ,t x )(kl + (x)u{t x , -oo)|0), (4-15) 

+ 0(i s/ -^)^|^_(y)«(t y ,^){r3^,-(x)r i n(t 1 .,-oo)|e), (4-17) 
d2 172 (x,y)=0(t x -g(0|{r3^,-(x)r i U (^,^){4 + (y)r 3 r 2 U (t y ,-oo)|0) 
+ 0(tj / -t a; )(0|{4 + (y)r3r 2 U (tj / ,t :c ){T3^,_(x)r i U (t :c ,-oo)|0). (4-18) 

Performing the thermal Bogoliubov transformations for Eq. (|4-4|) , we rewrite the 
thermal counter term Q by the transformed operators, ^ and £, 

Q = ~*J -0fMt*)&l- (4-i9) 

The thermal counter term satisfies = 0. This condition is important, and 

necessary to use the Feynman diagram procedure.® 

The thermal counter term can be fixed by the renormalization condition. Since 
Eq. (|4-19|) is proportional to the thermal counter term appears in the inverse 
propagator or the self energy, and modifies, at the leading order for the propagator, 
only d\ 2 and df in Eqs. g!5l) - (l4T8l ). Substituting Eq. gT9l) into Eq. (HT^ . we 
obtain 

<0|T[C(z)$(y)u(oo,-oo)]|0) - (e\T[fi(x)%(y)]\9) 
d 3 p 1 _ e -iu r (t x -ty) jp{x-y) 



(2tt) 3 2u p 

( fty ft* N 

x Q(t x - t y ) / dt s h p (t s ) + 9(t y - t x ) / dt s h p (t s ) 

\ J— oo J — oo 
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xB-\n p (t x ))( ° Q J )B(n p (t y )) 

r d 3 p 1 iuy.fa-tv) -ip-(x-y) 
J (27r)3 2w p 

(rt y rt x \ 

6{t x - t y ) I dt s h p (t s ) + 9{t y - t x ) I dt s h p (t s ) 
J — OO J —CO J 

XT 3 B(n p (t x )) T ( J o ) B-\n p {t y )) T n . (4-20) 

On the other hand, the quantum correction at the leading order is written in 
2x2 matrix form by using the free propagator, Dq 13 , and the self-energy, JJ 1112 , 

J d 4 z 1 d 4 z 2 D^ 1 (t x ,t zl ,x- Zl )iE^(t zl ,t Z2 , Zl - z 2 )D^(t Z2 ,t y ,z 2 - y) 

= J d 4 z 1 d 4 z 2 B' 1 (n lVxl (t x )) 

x f d I]r( x ~ zi)iU R {t Zl ,t Z2 ,z 1 - z 2 )D\} R {z 2 - y) SE B ,i(x,z 1 ,Z2,y) \ 
X V D o%( x ~ zi)i£A(t Z i,tz2,zi - z 2 )D 22 R {z 2 -y) ) 

xB(n^ yl (t y )) 

+ J cfV^B-Vlv.l^)) 

/ SS Bj2 (x,z 1 ,z 2 ,y) - D^ R (x - z x )iE R {t zl ,t Z2 ,z x - z 2 )D% 2 A (z 2 - y) 
\ -D 22 R (x - z 1 )%E A {t zl ,t Z2 ,z 1 - z 2 )Dl] A {z 2 -y) 

xB-\n^ y][ {t y )) T Tz 

+ J d A Zl d 4 z 2 T Z B{nyt x \{t x )) T 

x f - Dl] A {x - z 1 )iZ A (t Zl ,t Z2 ,z 1 - z 2 )D 22 R (z 2 - y) 

V - d o 2 a( x ~ ^^R^zi^zz^!- z 2 )D^ R (z 2 -y) 6E B ,3(x,z 1 ,z 2 ,y) 
xB(n^ yl (t y )) 

+ J d 4 Zl d 4 z 2 T 3 B(n^ xl (t x )) T 

x ( d I]a ( x - z i ) iE A {t zl ,t Z2 ,Z! -z 2 )D\} A {z 2 -y) \ 
\ 5E BA (x,z 1 ,z 2 ,y) Df A (x - zi)iE R {t zl ,t Z2 ,Zx - z 2 )Dl 2 A {z 2 - y) ) 

xB-\n^^t y )) T r^ (4-21) 

with 

5U B ,i(x,zi,z 2 ,y) 

= d I]r( x ~ 2i){ii7 12 (t Zl ,t Z2 ,zi - z 2 ) +iU R (t Zl ,t Z2 ,z 1 - z 2 )n^^(t Z2 ) 
-n lVzil (t zl )iU A (t Zl ,t Z2 ,z 1 - z 2 )}Dl 2 R {z 2 - y), (4-22) 
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5E B ,2(x,Z 1 ,Z2,y) 

= D^ R (x - Zi){ii: n (t Zl ,t Z2 ,zi - z 2 ) + iU R (t Zl - t Z2 ,z 1 - z 2 )n^^(t Z2 ) 

)}Dl] A {z 2 -y), (4-23) 

5U B , 3 (x,Z 1 ,Z 2 ,y) 

= ~ d o 2 a( x - z 1 )^iS 22 (t zl ,t Z2 ,z 1 - z 2 ) + iE R (t zl ,t Z2 ,z 1 - z 2 )n^^(t Z2 ) 

)}D 2 %(z 2 -y), (4-24) 

5S BA (x,zi,z 2 ,y) 

= Dl 2 A (x - ziJ^-iS 21 ^,^^! - z 2 ) - i£ R (t Zx ,t Z2 ,z x - z 2 )n^^(t Z2 ) 

+n\v zl \(t zl )iE A (t Zl ,t Z2 ,z 1 - z 2 )^Dl] A {z 2 -y). (4-25) 

An explicit form for the free propagator is given in (|A-3p . The retarded and advanced 
parts of the self-energy, S R and S A , are defined by 

S R = S 11 + S 12 = S 21 + S 22 , S A = U u - U 21 = U 22 - U 12 . (4-26) 

The first and last terms in the right-hand side of Eq. (|4-2ip have the same Bogoliubov 
transformation structure as the first and last ones in the right-hand side of Eq. ()4-20p , 
respectively. Below we identify Eqs. (|4-2(J|) and (|4-21|) as the thermal counter terms 
and the contribution of quantum corrections, respectively. 

H. Chu and H. Umezawa have proposed the self-consistency renormalization con- 
dition to fix the thermal counter term. The condition imposes {O\(,H,p(tx)£,H,p(t x )\0) = 
at the equal time limit, t x — > t y ^ where the subscript H denotes the Heisenberg 
picture and whose implication will be seen at the top of the next section. The condi- 
tion amounts to the vanishing off-diagonal elements, d\ 2 and d| , in the limit. Due 
to the tilde conjugation rules, both the equations give an equivalent condition. Thus 
the self-consistency renormalization conditions reduce to a single equation. From 
Eqs. (|4^20|) and Eq. (^2T|) we obtain 

d 3 p 1 



/*" dts [ -^—h j ts yp< x -y) 



+ lim [ d 4 z 1 d 4 z 2 5Z B ,i(x,z 1 ,z 2 ,y) = 0. (4-27) 

In thermal equilibrium 6Ub,i and 5X! B ^ vanish, so that these conditions are satis- 
fied automatically. In NETFD Eq. (|4-27|) shows the time evolution for the thermal 
Bogoliubov parameter, n p (t). 



Substituting the fields (|4-10[) - (14-13[) an d taking the equal time limit, we obtain 

lim df{x,y)=- f fP f* a _L J- e <p - iC e-'*-»(g|T[6 > (t,)fo(t.Hoo,-oo)]lg), 
•»->*» J (27r) J (27r) J y 2cj p y 2tjfe 

lim df(x,y)=- [ f* f* a * J- e-^ a5 e < *-»(g|T[6 > (t,)fo(t.Hoo,-oo)]lg). 
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For a practical calculation in homogeneous NETFD it is more convenient to 
employ the i-representation, while the spatial Fourier transformation is performed. 
Differentiating Eq. (|4-27|) with respect to t x and performing the spatial Fourier trans- 
formation, we obtain 

h p (t x ) = -2u p d t A lim dt Zl dt Z2 5HB,i(t x ,t Zl ,t Z2 ,ty;p)\, (4-28) 

where 6UB,i(t x ,t Zl ,t Z2 ,t y ;p) is written by the propagator and the self-energy in the 
t-representation, 



6U B ,i(t x ,t Zl ,t Z2 ,t y ;p) = Dl] R {t x - t Zl ;p){i£ 12 (t Zl ,t z . 2 ,p) 

P)}Dl 2 R (t Z2 -t y ;p). (4-29) 

This equation corresponds to the quantum Boltzmann equation for a relativistic 
neutral scalar field. 

§5. Boltzmann equation for a neutral scalar field 

The Heisenberg number density, n,H, p (t), is defined by 

(2^) 3 ^ 3 )(p - k)n H>p {t) = {e\a ] Hp (t)a Htk (t)\e). (5-1) 

It is the self-consistency condition®'^ in the previous section that establishes 
the correspondence between the thermal Bogoliubov parameter, n p (t), and the above 
particle number density, riH,p(t). We can confirm the correspondence as follows. 
The interaction hat-Hamiltonian, Hi, annihilates the bra vacuum (but not the ket 
vacuum in general though), 

{e\H! = 0, Hj\6) ± 0, (5-2) 

according to the thermal state conditions Eqs. (|2-18|) - ([2-2ip and the specific form of 
the thermal counter term Eq. (|4-19p . The ^-operators in the Heisenberg picture are 
defined by 

e H , p {t)=B(n p {t))^a% p {t), (5-3) 
% p it) = -4, p {t)B-\n p (t)f\ (5-4) 
with the Bogoliubov parameter n p (t) (not riH,p(t)) and satisfy 

= <^~L(*) = 0» W*)|0> + 0, Wt)|0> + 0. (5-5) 
Then it follows that 

(9\a^ p (t)a Hik (t)\9) 
= (0\ { W*) + (1 + np(t)) + n fc (t)4 jfc (t)} |0> 

= (e\£ H j,(t)tH,k(t)\o) + MtMUpmkM 9 )' ( 5 - 6 ) 
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X 



/ 



p \ , p 

< — 

k 2 

Fig. 1. 1-loop thermal self-energy derived three point scalar interaction model 



which implies 

(27rf5^(p - k)(n HtP (t) - n p {t)) = (0\£ Ba ,{t)£ H , k (t)\e) . (5-7) 

Thus the thermal Bogoliubov parameter is equal to the Heisenberg number distri- 
bution at each instant of time under the self-consistency renormalization condition, 
(0|£ff p(t)^H p(t)\0) = 0. In what follows we apply the self-consistency renormaliza- 
tion condition to the self-interacting systems of relativistic scalar field and derive the 
quantum Boltzmann equation for the thermal Bogoliubov parameter. 

5.1. A^ 3 interaction model 

First we calculate the time evolution of the thermal Bogoliubov parameter for 
a neutral scalar field with a three-point self-interaction. We start from the Hamilto- 
nian, 



H = H — H, (5-8) 



with 



H = d 6 x 



^KO) 2 + o (aO(-V* + m 2 )0 a (x)} + ^a(xf] , (5-9) 

where 4> a and 7r a mean cf>\ in Eq. (|3-15p and 7r„ in Eq. (|3-17p . respectively. The fields, 
cj) a and 7r a , are also equivalent to <% in Eq. (|3-16p and ff* in Eq. (|3-18|) . The tilde 
conjugate Hamiltonian, H, is described by the fields, cj> a and fr a , which are equivalent 
to 4>a an d ^a) respectively We evaluate the one- loop thermal self-energy by using 
the Feynman rules in the thermal doublet notationPS'^} J n the ^-representation we 
assign the propagator (|A-3p to each internal line and 

A Q = A ( \ \ , (5-10) 

to each vertex. 

At the one-loop level the thermal self-energy is diagrammatically represented in 
Fig. [H and is calculated as 



iE]jjL. loop (t Zl , t Z2 ; p) 

1 f d 3 ki d 3 k 2 , ^ , 



2 
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X 



2 2 2 

EE 



d z ki d 3 k 2 1 
(2tt)3 (2tt) 3 4w fel a; fc2 



(27r) 3 <5^(p-fe 1 + fc 2 ) 



il=l 12=1 ' 
-t^) e i(E kl , tl +E k2 , t2 )(t zl -t Z2 ) 



fki,ii{t Z2 ) fk 2 ,i 2 (t Z2 ) fkx,ii\Pz2) fj?2 ,*2 (^22) 

/fel ,ii (^22 ) /fc2 ,12 (^2 ) ~fki,ii(tz 2 ) fk2 ,12 (^22 ) 

+6>(t 22 e _i(i?fe i' l i+ £; fe2.*2)(^i-^2) 



/fcl ,ii (^21 ) fk 2 ,i 2 (t z 



fki ,ii (t zi ) fk 2 ,i2 \Pz\ ) 



fk%,i\ {tz\) fk 2 ,%2 {pz\ ) fkx,i\ (tzi)fk 2 ,i2 (^21 ) 



where 



-E^r,l — ^,2 



f q ,l(t) =n q (t), f q>2 (t) = l + n q (t), 
f q ,l(t) = l+n q (t), f q ,2(t) = n q (t). 



(5-11) 



(5-12) 
(5-13) 
(5-14) 



Thus the off-diagonal elements, (|4-22p and (|4-25|) . are given by 
lim / dt zl dt Z2 8U B ,i(tx,tz 1 ,t Z2 ,t y ;p) 
= lim / dt zl dt Z2 dE B ,4(t x ,t Zl ,t Z2 ,t y ;p) 

d 3 ki d 3 k 2 1 



1 ^-^-A 2 



dt s 



EE 



00 Pj 1= ii 2 =i 



(27r) 3 (27r) 3 4u] p uj kl uj k , 



(27r) 3 <5 (3) (p-fci + fc 2 ) 



sin{(a;p + ia + E ki2 ){t x -t s )} 

K 

Up + E kl ^ x + E k2 ^ 2 
^n p (t s )f kltil (t s )fk 2! i 2 (t s ) - (1 +np(t s ))f kltil (t s )fk 2> i 2 (t s )^, 



(5-15) 



at the equal time limit. When the Bose distribution function is assumed for n p (t), 
Eq. ()5-15p vanishes. It shows that the Bose distribution is a stationary solution for 
Eq. g2HD. 

Inserting Eq. (|5- 15[) into Eq. (|4-28p . we obtain the time evolution equation for 
the thermal Bogoliubov parameter, 



fl p (txj 



X 



tEE 



dt. 



d 3 fci d 3 k 2 



1 



(2tt) 3 (2vr) 3 4u p uj kl uj k2 



ii=l i 2 =l ' 

cos{(w p + E kuh + E k2ti2 ){t x - t,)}(2vr) 3 ^ 3 )(p - k x + fc 2 ) 

(%)(*s)/fci,ii(*s)/fc 2 ,ia(*s) - (! +»p(*s))/fci,ii(<a)/fe2,i 2 (*s)}- ( 5 ' 16 ) 
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The third line of this equation has the same statistical structure as the quantum 
Boltzmann equation. Off-shell contribution is included through the second line of 
Eq. QST5D . 

At the limit, t x — > oo, the second line of Eq. (|5T5p reduces to the delta function 
which guarantees the energy conservation. Then the collision term of the quantum 
Boltzmann equation is derived from Eq. (|5-15p . It shows that time evolution for the 
thermal Bogoliubov parameter is described by the quantum Boltzmann equation at 
the limit. 

5.2. A</> 4 interaction model 

Next we consider a neutral scalar field with a four-point self-interaction. The 
model is defined by the Hamiltonian, 

H = H — H, (5-17) 

with 

H = J d'x^TTaixf + M^i-Vl + m^M^ + ^M^]. (5-18) 

The numerical factor in the interaction term provides the same assignment to each 
vertex as in the previous model. 

We calculate the thermal self-energy in this model. Since there is no momentum 
transfer from the external to the internal lines, the self-energy has a diagonal form 
at the one-loop level. Thus the self-consistency renormalization condition is satisfied 
for hp{t) = at one-loop level. The time evolution of the thermal Bogoliubov 
parameter is induced from the two-loop self-energy illustrated in Fig. [2J We compute 
the diagram and obtain 

i^iij-ioopitz! , t Z2 ; p) 

1 f c2 3 fc! d 3 k 2 d 3 k 



3! 

T? S \ s {D (t Zl ,t Z2 ; fci)r3}^{^o(^,^ 2 ;fc2)r 3 } <572 {^o(^ 1 ,^ 2 ;fe3)r3} <572 A^ 
A 2 A A f d 3 k± d 3 k 2 d 3 k 3 1 3 (3) 



— ► 

/ N 
■> / ► \ ► 



P \ fe 2 J'P 



Fig. 2. 2-loop thermal self-energy in \(f> 4 interaction model 
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0(f z — t z ^ e i ( E k 1 ,n+ E k 2 ,i2+ E k 3 ,z :j )(iz 1 -t Z2 ) 

fk-L ,ll (i«2 ) fk-2 ,i2 (^22 ) fk 3 ,i 3 (tz 2 ) ~~ fkl ,i 1 (^22 ) f>!2 ,12 (^22 ) /fc3 ,13 (^22 . 

fki ,h (^2 ) /fca ,12 (^22 ) fk3,i3 (^2 2 ) ~~ /fci ,ii (^22) /&2 ,ia (^22 ) /fc3,«3 (^22 . 
/ fki,h\tzi)fk2, A3, 43(^21) — /fci,ji(^2i)/fc2, 12(^21 )/fc3,i 3 (i2i. 

V fkl,h\tzi) fk2 ,12 fyzi ) A3 , «3 (^Zl ) ~~ /fci ,il (^Zl ) /fc 2 ,12 (^21 ) /fe 3 ,13 (^2i . 

The off-diagonal elements, (|4-22p and (|4-25p . are given by 



(5-19) 



lim / dt Zl dt Z2 SE B ,i(t x ,t Zl ,t Z2 ,t y ;p) 

X^ty J 

= lim / dt Zl dt Z2 6E B ^(t x ,t Zl ,t Z2 ,t y ;p) 

tx—tty J 



** , 1 A 2 A A A /" d 3 fc! d 3 fc 2 d 3 k 3 1 



2w p 31 ( 27r ) 3 ( 27r ) 3 ( 27r ) 3 8^ fcl a; fc2 a; fc 3 



x sin(H + E kutl + E k2 , 2 + - t,)} (27r)3j( 3) (p _ fe _ fc2 _ fcj|) 

+ E kl ^ 1 + E k2 ^ 2 + E k3 ^ 3 

x np(ts)fki,h(t s )fk 2 ,i 2 (t s )fk 3> i 3 (t s ) - (1 +np(t s ))fk ltil (t s )fk 2 ,i 2 (t s )fk s> i 3 (t s ) 

(5-20) 

at the equal time limit. Substituting Eq. (|5-20|) to Eq. (|4-28p . we obtain the time 
evolution equation for the thermal Bogoliubov parameter, 

A 2 2 2 2 t d 3 fci d 3 fc2 d 3 fcs j 



A 2 /" 1 /" 

%>(y = c-i) ¥ E E E / dt W 



11=1 12=1 «3=1 



(27T) 3 (27r) 3 (27r) 3 %UJ p UJ kl UJ k2 U k3 

xcos{(w p + # fcliil + £ fc2) ; 2 + E k3)i3 ){t x - t s )}(2ir) 3 5 {3) (p - fci - fc 2 - fc 3 ) 

X Wp(t a ) fki,h(ts) fk2,i 2 ( t s) fk 3 ,i 3 (t s ) - (1 + n p (t s ))f kuil (t s )f k2: i 2 (t s )f kSi i 3 (t s ] 

(5-21) 

This equation has the same statistical structure as the quantum Boltzmann equation 
for the A(/> 4 interaction model. It should be noticed that a coefficient of the right-hand 
side in Eq. (|5*21|) is twice of the one obtained in Ref. [26]) . As is shown in Appendix 
B, Eq. (|5-21|) coincides with the quantum transport equation in the non-relativistic 
regime.^ 



§6. Conclusion 



We have investigated a relativistic neutral scalar field in NETFD. Thermal de- 
gree of freedom is introduced through the time dependent Bogoliubov transforma- 
tion. Then the thermal counter term has to be introduced for a consistent descrip- 
tion of both a(t) and £(t) under the transformation, and is a part of the interaction 
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Hamiltonian. We adopt the perturbative expansion, and calculate the full propaga- 
tor in the canonical formalism. The Bogoliubov matrix structure of the propagator 
is crucial. 

Applying the self-consistency renormalization condition 24 ' to a relativistic neu- 
tral scalar field, we have derived the time evolution equation for the thermal Bogoli- 
ubov parameter which is considered to be the particle number density. It has been 
shown that the equation reduces to the quantum Boltzmann equation in </> 3 and (j) 
interaction models. 

In this paper we impose the Lorentz covariance for the neutral scalar field and 
decompose it in terms of the creation and annihilation operators, £ and £>. It is 
not always possible to do so in a general situation of non-equilibrium system. Some 
modification would be necessary to apply the procedure to a field with a time- 
dependent screening mass, for example. 

There are some remaining problems. There is no counter term for the second 
and third lines in Eq. (|4T4p . These terms may have a nontrivial contribution to 
the time evolution equation at higher order.® In the present paper we have assumed 
spatial homogeneity. The space-time dependence is also important to study some 
relativistic systems. Some works to extend NETFD to spatially inhomogeneous 
systems have been attempted for non-relativistic field. An essence of such extension 
is to expand the field operator not by a complete set of plane wave functions, but 
by a complete set mixing momentum for diffusion process^ >EI3 or by a complete set 
of wave functions under trapping potential for cold atomic system while the 
equal-time commutation relations are preserved. We can formulate inhomogeneous 
TFD for relativistic fields in similar ways. We are also interested in applying the 
procedure to a relativistic Dirac field and an inhomogeneous system. We hope to 
solve these problems and report the result in future. 
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Appendix A 

Propagator for a free neutral scalar field 

In TFD the Feynman propagator for a free neutral scalar field is given by the 
expectation value of the time ordered product of two scalar fields. It has the 2x2 
matrix form in the thermal doublet notation,™ 

Df(t x ,t y ,x-y) ee (9\T[^(x)^ a (y)]\9), (A-l) 

where the neutral scalar fields, cj) a and 4> a , are decomposed into the positive and 
negative frequency parts by Eqs. (|3T5p and (|3T6p . The thermal Bogoliubov trans- 
formation is applied, then the propagator (|A-ip reads 
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+ {r 3 i?(n| V ^(^)) T }^M^C^)^ + (y)^){r3S- 1 (n |Vy| (^)fP^ 



+ {r 3 B(n lVxl (t x )Y} 



T\ 071 



72 



^172/3 



( y )^ 1 „(y)|0){r 3J B- 1 (n |Vy| (t 3/ )) T } 



(A-2) 



where the fields, and </>£,±, are defined in Eq. (|4-13p . 

Due to the definition of the thermal vacuum (|2-16|) and (|2-17|) the thermal prop- 



agator, Dq^ , reduces to 



(t x ,ty,x y) 



B-\n^{t x ))^DlJ(x-y)B(n^ Xt v )) 



72/3 



+{r,B{n\ Vx \{t x )) T r^DlJ(x - y){B-\n^ Xt y )) T n} 



72 



(A-3) 



where -Dj 1 ^ 2 (x — y) and -D^A 2 ( x ~~ 3/) represent retarded and advanced parts of the 
propagator, respectively, 



Dl] R {x-y) 
D 2 2 R (x-y) 
Dl] A {x-y) 
Dl 2 A {x-y) 



d 3 p 



6{tx — ty 



-ip-(x-y) 



(2vr)S • {y x) 2lo p 

—o{t y -t x )^-^-y\ 

V 



( 27r )3-v-« --2a; 



V 



(A-4) 
(A5) 
(A6) 
(A7) 



other components = 0. 



Thus we reproduce the result obtained in Refs. 2} and [22]) . The thermal propagator 
(|A 3|) has the same form as the one in an equilibrium system. The time dependence 
is introduced through the thermal Bogoliubov transformation. 



Appendix B 

Non-relativistic limit of the Boltzmann equation 

Here we take the non-relativistic limit of the time evolution equation for the A^ 4 
interaction model (|5-21|) and compare it with the transport equation for the cold 
atom system. 

The cold atom system is described by the Hamiltonian, 



H 



d 3 x 



V 2 g 
2m 2 



(B-l) 
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where ip represents a non-relativistic scalar field, 

The quantum transport equation for this system is given by^J 

d 3 p l d 3 p 2 d 3 p 3 



h p (t) = Ag 2 Re f ds f 

J — OO J 



(2tt) 3 (2vr) 3 (2vr) 3 
x(2tt) 3 5 {3 \ Pi +p 2 -p 3 -p) e - i ( e '>i +e ra-eP3-e P )(*-«) 

x{n pi (s)n P2 (s)(l +n P3 (s))(l + n p (s)) - (1 + n pi (s))(l + n p2 (s))n p3 (s)ra p (s)}, 

(B-3) 

where e p is the kinetic energy for the non-relativistic field, 

,2 



e p = — . (B-4) 



At the non-relativistic limit, \p\ <C m, the energy eigenvalue uj p reduces to 

ojp » m + e p . (B-5) 

We restrict the momentum for the scalar field in — p e < pi < p e with a cut-off 
parameter, < p e <C m. Thus the relativistic scalar field, <j> a , is decomposed to be 

M*)~ f ■^-^{a p (t x yP- x + ai(t x )e-*P- x } 

= a ,+ (x) + a ,_(x), (B-6) 

where 4> a:+ and <t> a - indicate the positive and negative frequency parts. 

The interaction for the cold atom system (jB- lj) can be identified with the inter- 
action, (p at +<p a}+ (p a ^(f) a ^. There are six corresponding terms in the interaction, A(%. 
Thus we obtain the following correspondence between the interaction terms for the 
cold atom system and the Xcfr* model in the non-relativistic regime, 

^ (x)^ (x)i){x)i){x) =6^0 a)+ (a:)0 a)+ (x)<^ o _(x)0 o _(a;), (B-7) 

We find that there is a correspondence if we make a replacement 

A O 8m 2 g. (B-8) 

Since the transport equation ()B-3p comes from the two-body scattering, we pick up 
terms which represent the two-body scattering in Eq. (|5-21|) . Hence Eq. (|5-21|) is 
rewritten as 

A 2 /■'» , f d 3 ki d 3 k 2 d 3 k 3 1 



n p (t x ) = — J dt s J 



(2tt) 3 (2tt) 3 (27r) 3 8u p uj kl uj k2 uj k 
xcos{(o; p + u kl - oj k2 -u) k3 ){t x - t s )}(2vr) 3 (5 (3) (p- ki - k 2 - fe 3 ) 

x n p (t s )n kl (t s )(l + n k2 (t s ))(l + n ks (t s )) - (1 + n p (t s ))(l + n kl (t s ))n k2 (t s )n k3 {t s ) 

(B-9) 
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Substituting (|B-5P and (|B-8|) into Eq. ()B-9[) we obtain the quantum Boltzmann equa- 
tion in the non-relativistic regime, 



4<r 



I dt s f 

J -co J 



d 3 k t d 3 k 2 d 3 k 3 



(2tt) 3 (2vr) 3 (2tt) 3 
xcos{(e p + e fcl - e k2 - e k , A ){t x - i s )}(2vr) 3 J (3) (p - k x - k 2 - k 3 ) 

n p (t s )n kl (t s )(l + n k2 (t s ))(l + n ks (t s )) - (1 + n p (i s ))(l + n kl (t s ))n k2 (t s )n k3 (t s ) 

(B40) 



In homogeneous system the quantum Boltzmann equation (|B-10j) and Eq. ()B-3j) 
become identical. Therefore the time evolution equation (|5-2ip is consistent with 
the quantum transport equation (|B-3p . 
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